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We present lattice studies of heavy quark potentials in the quenched approximation of QCD at finite tempera- 
tures. Both, the color singlet and color averaged potentials are calculated. While the potentials are well known at 
large distances, we give a detailed analysis of their short distance behavior (from O.Of 5 fm to f fm) near the critical 
temperature. At these distances we expect that the T-dependent potentials go over into the zero temperature 
potential. Indeed, we find evidences that the temperature influence gets suppressed and the potentials starts to 
become a unique function of the underlying distance scale. We use this feature to normalize the heavy quark 
potentials at short distances and extract the free energy of the quark system in a gluonic heat bath. 



1. INTRODUCTION 

The study of heavy quark potentials (free en- 
ergies) is important for our understanding of the 
fundamental forces in strongly interacting mat- 
ter. For example it finds applications in the 
phenomenology of heavy quarkonium physics at 
finite temperature [1,2,3]. Our aim is to get 
deeper insights in their short distance behavior 
{R < 0.1/m). It is expected that at distances 
R << 1/T the presence of a thermal medium 
does not modify physics. However the question 
arises at which distance scales the temperature 
effects become negligible. 

We performed lattice studies of the color aver- 
aged and the color singlet heavy quark free en- 
ergy. Our calculations have been performed on 
a 32'^ X 8 and 32'^ x 16 lattice using a tree-level 
Symanzik improved action and on a 64^ X 16 lat- 
tice using the standard Wilson action of the pure 
S't/(3) theory. Up to an additive renormalization, 
the color averaged free energy, Fav '■= Fav{R,T), 
is defined in terms of Polyakov loop correlation 
functions [4] Fay = -Tln(frL(0)frLt(i?)) +(7^^, 
where L denotes the Polyakov loop and Tr is 
the normalized trace: Tr = 1/A^cTr. Follow- 
ing Ref. [5] the color singlet free energy Fg :— 
Fg {R, T) is identified with the cyclic Wilson loop 
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via Fs = -T \n{f TL{0)Dj^L^R)Dl) + C,, where 
Dj^ is the product of links on the shortest con- 
nection from to R. This formula is expected to 
be valid at short distances [5]. The renormaliza- 
tion constants Cav and Cg contain temperature 
independent divergent self-energy contributions. 
As no additional divergences are introduced at 
T > we can fix the renormalization constants, 
if we are able to separate the T-independent part 
in the free energy and match it to the zero tem- 
perature heavy quark potential. 

The color averaged free energy can be related 
to the singlet (Fg) and octet {Fg) contributions 
due to: 

exp(-F,,/r) 

- 1/9 exp(-F,/T) + 8/9 exp(-Fo/r). (1) 

In a static system the internal energy reduces to 
the potential V and thus the free energy refers to 
F — V — TS, where S denotes the entropy. In 
a perturbative treatment of QCD at lowest order 
the singlet and octet potentials are of Coulomb 
form [6]: 

ysiR)--l£i^ VoiR)^~lvAR), (2) 

where g is the QCD running coupling. A high 
temperature expansion of (1) in combination with 
(2) yields the color averaged heavy quark poten- 
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tial: 

V{R,T) ^ _1 l_ ^ al 

T 9 (47r)2 {RTf %RTf ' ^ ' 

with as = g'^ /{Ait). Note that (3) is restricted to 
large distances in order to justify the high temper- 
ature expansion. At short distances the octet con- 
tribution in (1) is suppressed leading to a singlet- 
like behavior of the color averaged potential [7] . 

2. THE COLOR AVERAGED POTEN- 
TIAL AT T > 

We study the short distance behavior of the 
color averaged potential in terms of the screen- 
ing function S{R,T) defined by S{R,T) = 
-W{R,T)T[1/R]-'^, with being the lattice 
Coulomb potential. At large distances the screen- 
ing function is exponentially suppressed (hence its 
name) and in the regime where the perturbative 
expression (3) is applicable we expect S = a^. 
In Fig. 1 the square root of the screening func- 
tion is plotted. Two distinct features arc evident 
from our data: First, the screening function is 
decreasing at very short distances (i?r<0.1) indi- 
cating the singlet-like behavior mentioned before. 
Secondly, screening effects become important at 
rather short distances (about 0.3/T). Nonethe- 
less, there still remain temperature effects in the 
color averaged potentials at short distances. 

We have also compared the lattice data with 
a 1-loop perturbative calculation (up to 0{g^)) 
in Fig. 1. In this calculation we have added the 
non-static contributions to the potential which 
are expected to become important in the short 
distance regime (see Appendix A). As one can 
see from Fig. 1 the color averaged potential can 
be described by 1-loop perturbation theory for 
T>6Tc and distances 0.3<iiT<l. For smaller 
temperatures and larger distances perturbation 
theory breaks down due to non-perturbative ef- 
fects. Our calculation also breaks down at very 
small distances RT<Q.l because of the enclosed 
high temperature expansion. 




0.05 I • • • ■ • • ■ • 1 

0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 



RT 

Figure 1. The square root of the screening func- 
tion \/S calculated on a 32"' x 8 lattice (filled sym- 
bols) and on a 32^^ x 16 (or 64"^ x 16) lattice (open 
symbols) and the corresponding perturbative re- 
sults (lines). 



3. THE TEMPERATURE INDEPEN- 
DENCE AT SHORT DISTANCES 

In order to eliminate effects from color averag- 
ing we now refer to the short distance behavior of 
the color singlet potentials. Our results together 
with the zero temperature potential from [8] are 
shown in Fig 2. The potentials have been nor- 
malized at short distances to make a comparison 
with the zero temperature limit possible. As one 
can see from Fig. 2 the temperature dependence 
of the heavy quark potential becomes first visi- 
ble at distances of about R>0.1fm. Thus the 
(finite T) color singlet potentials approach their 
zero temperature limit already at the plotted dis- 
tances. Indeed this analysis provides evidence 
that the running coupling extracted from color 
singlet data using (2) starts to get T-independent 
at these distances. 

We now argue that these normalized data refer 
to the free energy of the estimated quark system - 
enclosing effects from the gluonic heat bath - be- 
cause divergent self-energy contributions cancel 
out. Moreover we expect that the free energy has 
a well defined continuum limit in our approach 
because the normalization gets independent of 
the lattice cut-off at sufficiently short distances 
(see Fig. 2). 
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Figure 2. The zero temperature limit of color 
singlet potentials at non-zero temperature. 



4. CONCLUSIONS 

Our study shows that the heavy quark poten- 
tials have a non-trivial structure at short dis- 
tances. Color averaging influences the physics of 
potentials at short distances and causes that the 
averaged data approach the T = limit at shorter 
distances as do the singlet data. The color singlet 
potentials approach the zero temperature limit at 
distances R ^ 0.1/m and T ^ Tc- Due to this 
feature we have extracted the free energy of a 
heavy quark pair enclosing the gluonic heat bath. 
A comparison of our 1-loop perturbative calcula- 
tion with lattice data suggests that the pertur- 
bative expansion is applicable for temperatures 
T>6Tc and distances 0.3<i?r<l. 



THE NON-STATIC 
TIONS 



CONTRIBU- 



We have calculated the color averaged heavy 
quark potential at finite temperature based on 
S[A,, Ao] = / dzo / <Px{-l/2{FS,f - 1/4(^^5)2). 
In this expression F^^ is the usual non-abelian 
field-strength tensor where Ai {i = 1,2,3) de- 
notes the static and Aq the non-static gauge field 
and we introduce a temporal gauge (i9oAo = 0). 
The (7^-correction (1-loop) to the color averaged 
potential can be divided in a static (A]/**) and a 
non-static (AV^"*) part: 



The static contribution is well known from [5]. 
For the non-static part we find at 1-loop level: 

Al^""(i?,T) = 
g4/32(]v2 - 1) /e-™o« 
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(5) 



where IIqq [p) is the self-energy of the non-static 
field with \\^^{p 0) to|) and /3 = l/T. Note 
that in our conventions the leading order term 
comes from the non-static contribution. In fi- 
nite temperature field theory II^q{p) consists of 
the UV-divergent zero temperature part and the 
convergent thermal part: 11'^^ {p) — IIqq {p) + 



nop (p) • IIqo {p) is calculated in [9] : 
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11 , 47rz/2e-T^ 
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(6) 



The thermal contribution can be extracted from a 
direct integration of the loop diagrams contribut- 
ing to I{^q{p) with standard calculation tech- 
niques at finite temperature. Our 1-loop calcu- 
lation of IIqq (p) is in accordance with (6) and the 
resulting potentials are plotted in Fig. 1. 
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^Vav{R,T) = AV''*{R,T) + AV'''{R,T). (A) 



